function [s_times, s_val] = pareto(max_time, onalpha, offalpha, onmu, offmu)

nproc = 1;      % number of processes to generate
on_par = {onalpha, 1/(onmu*(onalpha-1))};
off_par = {offalpha, 1/(offmu*(offalpha-1))};
[jmptimes, indproc]=onoff(nproc, max_time, @simparetonrm, on_par, @simparetonrm, off_par, 0);

[s_times, s_val] = stairsum(jmptimes, indproc);

zeri = find(s_times == 0);
s_times = s_times(zeri(end):end);
s_val = s_val(zeri(end):end);

zeri = find(s_times == max_time);
s_times = s_times(1:zeri(1));
s_val = s_val(1:zeri(1));

%clf;
%hold on;
%stairs(jmptimes, indproc, 'b');
%stairs(s_times, s_val, 'r');
%hold off;

end